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. Aims. We present the stellar mass functions (SMFs) of star-forming and quiescent galaxies from observations of 10 rich, red- 

Si^ I sequence selected, clusters in the Gemini Cluster Astrophysics Spectroscopic Survey (GCLASS) in the redshift range 0.86 < z < 

1.34. We compare our results with field measurements at similar redshifts using data from a Kj-band selected catalogue of the 
COSMOSAJltra VISTA field. 

Methods. We construct a Kj-band selected multi-colour catalogue for the clusters in 10 photometric bands covering M-8yum, and 
Q I estimate photometric redshifts and stellar masses using SED fitting techniques. To correct for interlopers in our cluster sample, we 

5— ( , use the deep spectroscopic component of GCLASS, which contains spectra for 1282 identified cluster and field galaxies taken with 

f/j ■ Gemini/GMOS. This allows us to correct cluster number counts from a photometric selection for false positive and false negative 

Ci ' identifications. Both the photometric and spectroscopic samples are sufficiently deep that we can probe the stellar mass function down 

to masses of 10'" Mq. 

Results. We distinguish between star-forming and quiescent galaxies using the rest-frame U-V versus V-J diagram, and find that the 
best-fitting Schechter parameters a and M* are similar within the uncertainties for these galaxy types within the different environ- 
^ ' ments. However, there is a significant difference in the shape and normalisation of the total stellar mass function between the clusters 

and the field sample. This difference in the total stellar mass function is primarily a reflection of the increased fraction of quiescent 
galaxies in high-density environments. We apply a simple quenching model that includes components of mass- and environment- 
^ft ■ driven quenching, and find that in this picture 44^j% of the star- forming galaxies, which normally would be forming stars in the field, 

l*^ ■ are quenched by the cluster 

^h ' Conclusions. If galaxies in clusters and the field quench their star formation via different mechanisms, these processes have to con- 

f^ , spire in such a way that the shapes of the quiescent and star-forming SMF remain similar in these different environments. 

" ] , Key words. Galaxies; clusters: general - Galaxies: mass function - Galaxies: evolution - Galaxies: photometry 

k>( 1 . Introduction els, and these models have to match and predict the SMF for a 

'^_j [ range of redshifts and environmental densities. 

.'^: One of the missing parts in the theory of galaxy formation Although models are tuned to match the observations at z=0, 
and evolution is a detailed understanding of the build up of there is in general still a poor agreement between observations 
stellar mass in the Universe. While the hierarchical growth of and theory at higher redshift. Models generally show an ex- 
dark matte r haloes has been st udied in large N-body Simula- cess of galaxies with a stellar mass (M*) ~ 10'° Mp around 
tions (e.g. LSoringel et al.ll2005l) . the baryonic physics that reg- z -1-2 compared to o bservational data (e.g. [Bower et al.||20ll 
ulates the cooling of gas and formation of stars in these haloes IWeinmannet alj|2012|). At higher redshifts the number of high- 
is much harder to simulate and is not yet well understood. To mass galaxies is generally underpredicted by the models. For a 
understand which physical processes are dominant in shaping detailed comparison between models and the observed SMF also 
the stellar content of galaxies, models need good observational see |Marchesini et al. ( 2009) . 

constraints. One of the most fundamental observables of a pop- At low redshifts (z < 0.2) the SMF has been measured 

ulation of galaxies is their stellar mass function (SMF), which from wide field data and spectroscopic information (ICole et al.l 

describes the number density of galaxies as a function of stel- 2001; B ell et al. 2003) , while at higher redshifts the SMF has 

lar mass. Measuring the SMF as a function of cosmic time pro- been measured fro m deep surveys by majking use of photometric 

vides useful constraints on the parameters in semi-analytic mod- redshift estimates dPerez-Gonzalez et al.ll2008t iMarchesini et al.l 
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l2009t lllbert et aDl2010l) . The general consensus is that the to- 
tal stellar mass density evolves slowly between < z < 1, 
which can also be inferred from the sharp decline of cosmic 
star formation in the Lilly-Madau diagram (iLillv et al.lll996t 
iMadau et al.ll 19961) in this redshift range. The main evolution is 
in the normalisation of the SMF, where as the shape does not 
show a substantial evolution since z ~ 4 dPerez-Gonzalez et al.l 
|2008|) . However, since these deep surveys generally probe small 
volumes, the do minant source of random uncertainty is often 
cosmic variance dSomerville et al.l l2004t IScoville et al.l l20()7t 
iMarchesini et al.ir2009l) . which is expected to not only have an 
eifect on the norma lisation but also on the shape of the observed 
SMF (iTrenti & Stiavelli 2008) . Observations over large areas, or 
a combination of multiple sight lines, are used to reduce this 
source of uncertainty. 

Besides the general time evolution of the properties of galax- 
ies, they are also observed to be strongly influenced by the den- 
sity of their environment. In particular, galaxies in overdense 
regions show lower star formation rates, and a higher fraction 
of red galaxies. At low redshifts, the Sloan Digital Sky Survey 
has allowed us to quantify these correlations with high precisio n 
dKaufi^mann et al.l2004tlBalogh et alJ2004tlBlanton et al.l2005h . 
The fraction of galaxies that are red is also a function of their 
stellar mass, with more massive galaxies being redder and form- 
ing fewer stars. The quenching fraction of galaxies being a func- 
tion of both stellar mass and environmental density, some recent 
studies have suggested the processes of "mass quenching" and 
"environmental quenchi ng" to be operating completely iridepen - 
dently from each other (iPeng et al.ll201ol iMuzzin etal.1 120121) . 
each operating on different time scales. Naively, we would ex- 
pect the combination of these processes to affect the shape of the 
SMF. 

A measurement of the SMF of galaxies as a function of en- 
vironmental density therefore provides further constraints on the 
physical processes that are important in these dense regions. For 
example, galaxies falling into massive galaxy clusters are ex- 
pected to be stripped of their cold gas component due to ram- 
pressure stripping, and a lack of new inflowing cold gas leads 
to a galaxy's star formation being turned off. Galaxies in groups 
and clusters are also expected to interact gravitationally through 
mergers and experience strong tidal forces as they fall towards 
the cluster centre. 

Combining these measurements done over a range of red- 
shifts and environments puts constraints on the way galaxies 
quench their star formation, since it allows one to separate be- 
tween internally and externally driven processes. Some stud- 
ies have attempted to measure the SMF a s a function of lo- 
cal environment at 0.4 < z < 1.2 (e.g . Bundv et al. '2006'; 
Bolzonella et al.ll2010t PVulcani et al.ll201 iL 2012: Giodini et al.i 



2012h . A measurement of the SMF at the highest densities has 



not yet been achieved in this redshift range. This is partly be- 
cause the deep (and therefore limited in area) surveys that have 
been used for SMF measurements (mostly the COSMOS and 
DEEP2 fields) do not contain the extreme overdensities corre- 
sponding to the most massive clusters of galaxies. 

In this paper we present a measurement of the SMF of galax- 
ies in 10 rich galaxy clusters at a range of redshifts (0.86 < 
z < 1.34). These clusters were detected using the red-sequence 
method on data from the Spitze r Adaptation of the Red-sequence 
Clust e r Survey (SpARCS, s ee iMuzzin et all2009HWilson et al.l 
l2009t iDemarco et al.l 120101) . and have typical velocity disper- 
sions of cr,, = 700 km/s which imply halo masses of M200 - 
3 X 10'"^ Mo. We combine deep photometric data in 10 or 11 
bands with the extensive deep spectroscopic coverage that we 



obtained from the Gemini Cluster Astr ophysics Spectroscopic 
Survey (GCLASS, IMuzzin et all 1201 2h . This allows us to es- 
timate stellar masses for individual objects and quantify the 
amount of interlopers in the photo-z selected sample as a func- 
tion of mass and projected clustercentric distance. We use the 
UVJ-diagram to photometrically separate between star-forming 
and quiescent galaxies, which is critical because the two galaxy 
types suffer from different observational difficulties and com- 
pletenesses. We also provide a comparison between the clus- 
ter results and the SMF measured from Ultra VISTA/COSMOS 
field. 

The structure of the paper is as follows. In Sect. |2]we give an 
overview of GCLASS, and the spectroscopic and photometric 
data that have been taken for this cluster sample. We also de- 
scribe the data from the reference UltraVISTA survey. In Sect.|3] 
we present our measurements of photometric redshifts, stel- 
lar masses and rest-frame colours to distinguish between star- 
forming and quiescent galaxies. We also explain how we correct 
the photometric sample for incompleteness by making use of the 
spectroscopic data. In Sect. |4] we present our results and make 
comparisons between the two galaxy types, and between clus- 
ter environments and the field. In Sect. |5]we discuss our results 
in the context of galaxy evolutionary processes and in particu- 
lar quenching in these massive clusters. We summarise and con- 
clude in Sect. |6] Extra information considering colour measure- 
ments and calibration are presented in the Appendices. There we 
also compare the UltraVISTA field SMF with the field probed by 
GCLASS outside the clusters to test for possible systematics. 

All magnitudes we quote are in the AB magnitudes system 
and we adopt ACDM cosmology with Oni = 0.3, Qa = 0.7 and 
Ho = 70kms-iMpc '. 



2. Sample & Data description 

2.1. The GCLASS cluster sample 

The GCLASS cluster sample consists of 10 of the richest clus- 
ters in the redshift range 0.86 < z < 1.34 selected from the 
42 square degree SpARCS survey, see Table |T] Clusters in the 
SpARCS survey were detected using the cluster red-s equence 
detection method developed by iGladders & Yed (l2000h . where 
the z' - 3.6yum colo ur was used to samp le the 4000A break at 
these redshifts (see IMuzzin et al.l |2008|) . For an ex tended de- 
scription of the SpA RCS survey we refer to Muzzin et al. (20091). 
IWilson et al] (|2009|) and iDemarco et al.1 (|2010|) . The 10 clusters 
that were select ed from the SpARC S survey for further study 
are described in IMuzzin et al.l (1201 2l) . and can be considered as 
a fair representation of IR-selected rich clusters within this red- 
shift range. 

We note that there is a possible selection bias in favour 
of systems with a high number of bright red galaxies. It is 
impossible to select clusters based on their total halo mass 
and therefore any cluster sample has potential selection bi- 
ases, whether it is X-ray selected, SZ-selected, or galaxy- 
selected. We note that follow-up studies of X-ray or SZ- 
selected clusters in the same redshift range also show a signifi- 
cant oyer^densitYofred-sequence galaxies (e.g. Blakeslee et al] 
I2OO3I: lMunisetal.1120051) Furthermore, the field SMF at z=l 
shows (e.g. IMuzzin et al.l l2013al) that even in the field, the 
bright/massive end of the population is completely dominated 
by red galaxies. Therefore it seems unlikely that a red-sequence 
selection results in a significant selection bias, at least for the 
most massive clusters at a given redshift such as the GCLASS 
sample. 
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Table 1. The 10 GCLASS clusters selected from SpARCS for follow-up spectroscopic and photometric observations. These clusters 
form the basis of this study. 



Name" 



RA 
J2000 



Dec 

J2000 



Ks-band IQ 
PSF FWHM ["] 



Klim 

[magAfl] 



[Mo] 



limit from bc03' 



SpARCS-0034 


0.867 


00:34:42.086 


-43:07:53.360 


1.01 


21.53 


10.42 


10.43 


SpARCS-0035 


1.335 


00:35:49.700 


-43:12:24.160 


0.40 


23.60 


9.92 


9.95 


SpARCS-0036= 


0.869 


00:36:45.039 


-44:10:49.911 


1.23(J) 


22. 11 (J) 


10.53 


10.50 


SpARCS-0215 


1.004 


02:15:23.200 


-03:43:34.482 


1.00 


21.73 


10.45 


10.46 


SpARCS-1047f 


0.956 


10:47:32.952 


57:41:24.340 


0.61 


22.68 


10.17 


10.04 


SpARCS- 105 If 


1.035 


10:51:05.560 


58:18:15.520 


0.86 


22.96 


9.99 


9.99 


SpARCS-1613 


0.871 


16:13:14.641 


56:49:29.504 


0.81 


22.55 


9.97 


10.02 


SpARCS-1616 


1.156 


16:16:41.232 


55:45:25.708 


0.84 


22.65 


10.33 


10.20 


SpARCS-1634 


1.177 


16:34:35.402 


40:21:51.588 


0.77 


22.88 


10.14 


10.13 


SpARCS-1638 


1.196 


16:38:51.625 


40:38:42.893 


0.66 


23.00 


10.13 


10.09 



For full names we refer to lMuzzin et al.l(l2012h . 
80% completeness limit for simulated sources. 

Corresponding mass completeness limit based on the galaxy in Ultra VISTA with the highest M/L fitted at thi s redshift at Kum. 

Mass limit from a synthetic spectrum with t = lOMyr starting at age of universe at that redshift with no dust dBruzual & Charlodl2003^ . 

For SpARCS-0036 we used to J-band as the selection band since it is signiiicantly deeper than the Ks-band. The image quality and magnitude 

limits refer to the J-band for this cluster 

Since the BCG is oifset from the centre, this is a better approximation for the cluster centre (different from lMuzzin et al.](l2012h ). 



Table 2, Properties of the 10 GCLASS clusters. 



Name 


^spec 


(Tv 


M200 


-^^200 






[km/s] 


[lO'^Mol 


[Mpc] 


SpARCS-0034 


0.867 


700!^" 
780!| 
750! 
640!;| 


3 6+' * 


1.1+"-' 
9+^-^ 

o.9!|^ 


SpARCS-0035 


1.335 


SpARCS-0036 


0.869 


45+1.6 
^•"-1.3 


SpARCS-0215 


1.004 


SpARCS-1047 


0.956 


(SfiO+'O 

1350!; 


9 Q+1.0 
J 2+6-3 

26.1+^1 
2.8!;"' 


1,0+01 

0.7+^^-^ 
2 i+8:i 

9+^-^ 
"■''-0.1 


SpARCS-1051 


1.035 


SpARCS-1613 


0.871 


SpARCS-1616 


1.156 


SpARCS-1634 


1.177 


7gr)+60 


44+11 

1 o+'''-^ 


1-011 


SpARCS-1638 


1.196 


480+50 

too_|„n 


0.6!- 



2.2. Spectroscopy 

The clusters in the GCLASS sample have extensive optical spec- 
troscopy, which has been taken using the GMOS instruments on 
Gemini-North and -South. For details on these measurements, 
the target select ion and an overview of the reduction of these 
data, we refer to lMuzzin et al. (2012) . 

In summary, spectroscopic targets were selected using a 
combination of their 3.6fim fluxes, z' - 3.6yL(m colours, and their 
projected clustercentric radii. The colour priority selection was 
chosen to be sufficiently broad so that there is no selection bias 
against blue galaxies within the cluster's redshift range. Because 
the mass-to-light ratio in the 3.6 jum channel is only a weak func- 
tion of galaxy type, the targeting completeness is, to first order, 
a function of radial distance and stellar mass only. The assigned 
targeting pr iority is highest for massive objects near the cluster 
centres (see lMuzzin et al.ll2012i Fig. 4). 

For these 10 clusters there are 1282 galaxies in total with 
redshifts, of which 457 are cluster members. For more than 90% 
of the targeted objects with stellar masses exceeding IO'^Mq, 
the limiting mass of the photometric data, a redshift was mea- 
sured with high confidence. Note that the targeting prioritization 



is known, we do not select against a particular type of galax- 
ies, and we have a high success rate of obtaining redshifts over 
the stellar mass range we study. Therefore, although the spec- 
troscopic sample is incomplete, it is a representative sample for 
the underlying population of cluster galaxies. The targeting com- 
pleteness can be quantified, and in Sect. 13.41 we use the spectro- 
scopic sub-sample to correct the full sample for cluster member- 
ship. 

We have performed a dynamical analysis (Wilson et al., in 
prep) to study the distribution of line-of-sight (LOS) velocities 
of the spectroscopic targets. For all 10 clusters, the distribution 
of LOS velocities approximates a Gaussian profile, which is an 
indication that the clusters are (close to) virialised. From this 
distribution we measure the LOS velocity dispersion (o",,) of the 
clusters. This leads to estimates of /?200, the radius at which the 
mean interior density is 200 times the critical density (pcnt), and 
M200, the mass contained with i n j^200 - The current analysis is 
done similar to iDemarco et al.l ([2010), and is based on an ex- 
panded spectroscopic data set. Table|2]shows the cluster proper- 
ties obtained from this analysis. 

The clusters have typical velocity dispersions of cr,, - 700 
km/s which imply halo masses of M200 - 3 x IO'^^Mq. Note that 
SpARCS-1613 is much more massive, with a velocity dispersion 
of cr,, = 1350 km/s. This high value is consistent with the X-ray 
temperature measured from a recent Chandra observation (see 
ElUngson, in prep.). 

2.3. Photometric Data 

Optical ugriz data for the six clusters observable from the 
Northern sky were taken with MegaCam at the Canada-France- 
Hawaii Telescope (CFHT). For the clusters in the South, gri 
data were taken with IMACS at the Magellan telescopes, 
and the z-band data using the MOSAIC-II camera mounted 
on the Blanco telescope at the Cerro Tololo Inter-American 
Observatory (CTIO). There is J- and Kj-band imaging data from 
WIRCam at the CFHT for the Northern clusters, and from ISPI 
at the Blanco telescope or HAWK-I at the Very Large Telescope 
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(VLT) UT4 for the Southern clusters. Note that these near- IR 
data were already presented and used in iLidman et al.l (1201 2l) to 
study the evolution of BCGs. The photometric data set also in- 
cludes the 3.6, 4.5, 5 8 and S.Oyum IRAC channels from SWIRE 
jLonsdale et al.l2003l) . For half of the clusters, including the four 
at the highest redshifts, we obtained deeper IRAC observations 
from the GTO programs PID:40033 and PID:50161. The mea- 
sured depths and an overview of instruments that were used are 
listed in the Appendix in Table I A. II Note that for three of the 
Southern clusters we do not have M-band data coverage, but this 
typically only affects photometric redshift estimates for z < 0.3. 
In Appendix lAl we give a comprehensive description of the 
photometric data processing leading to a multi-wavelength cov- 
erage with a field of view of 10' x 10' centred on the Southern 
clusters, and a 15' x 15' field of view for the Northern clusters. 
This wide field view provides information up to several Mpc 
from the cluster centres at the respective cluster redshifts, even 
for clusters at the lowest redshifts. 



2.3.1. Object detection 

To measure the stellar mass function it is necessary to obtain 
a catalogue in which the galaxy sample is complete down to a 
known mass threshold, independent of their star- formation prop- 
erties. In an IR-wavelength band the M/L varies little for differ- 
ent star formation histories, so that the luminosity in those bands 
is a good tracer for the total stellar mass of a galaxy. 

Because the IRAC channels suffer from a large point spread 
function (PSF), separating objects on the sky is difficult. As a 
compromise between good image quality and detection in a red 
filter, we therefore choose the Kj-band as the selection band. 
We use SExtractor to detect all sources in the Ks-band that have 
5 adjacent pixels with significance > 2.5cr relative to the local 
background. 

We obtain a clean catalogue by excluding regions near bright 
stars and their diffraction spikes, and separate stars from galax- 
ies by using their observed colours. In the u - J versus J - K 
colour plane the distinction between stars and galaxies is clear 
(see e.g. Whitakeret al. 2011), and we find that the following 
colour criterion can be used to select a sample of galaxies. 

J - K > 0.18 ■ (m - J) - 0.70 U J - K > 0.08 ■ (u - J) - 0.40 (1) 

When M-band data are not available (as is the case for 3 clus- 
ters), we use the g - i colour to select galaxies by 

J - K > 0.19 ■ (§ - J) - 0.45 U J - K > 0.07 ■ (g - J) - 0.25. (2) 



2.3.2. Colour measurements 

To measure photometric redshifts and stellar masses for the 
galaxies, accurate colour measurements are necessary. The cir- 
cumstances of the atmosphere and optical instruments change 
continuously, and therefore the shape and size of the PSF is dif- 
ferent between telescopes, exposures and filters. Therefore it is 
non-trivial to measure colours of the same intrinsic part of a 
galaxy. A common approach is to degrade the PSF of the im- 
ages in all filters to the PSF of the worst seeing, after which the 
colours are measured by comparing the flux in fixed apertures 
for all filters. 

An alternative approach, proposed by iKuiikenI (l2008h . is 
to perform a convolution of the images in each filter with 
a position-dependent convolution kernel to make the PSF 
Gaussian, circular and uniform on each image. The images in the 
different filters are not required to share the same Gaussian PSF 



size. Fluxes are measured in apertures with a circular Gaussian 
weighting function, whose size is adapted to ensure that the same 
part of the source is measured. Because the weighting function 
approximately matches the galaxy profiles, this technique gen- 
erally improves the S/N of the measurement compared to a nor- 
mal top-hat shaped aperture, and we elaborate on this method 
in Appendix lA.2l Note that this technique is not suited for mea- 
surements of the total flux, only to obtain colours of a galaxy. 

The photometric zeropoints are set based on standard-star 
observations. We improve the precision of the zeropoints in the 
ugrizJKfi filters by iiiaking use of the universality of the stel- 
lar locus dHigh et al.ll2009l) and comparing the mea sured stellar 
colours in our images with a reference catalogue (ICovev et al.l 
|2007|) . Further details can be found in Appendix lA. 21 

2.4. UltaVISTA field reference 

In this paper we compare the cluster results to measurements 
from a representative field at similar redshifts as the clusters. 
The last decade has seen substantial improvement in the depth 
and an increase in the field-of-view of ground-based NIR sur- 
veys. The most recent wide-field NIR survey is UltraVISTA 
dMcCracken et al.l2012h . which is composed of deep YJHKj im- 
ages taken using the VISTA telescope on a 1.6 square degree 
field that overlaps with COSMOS. 

The field sample in this study is based on a Kj- 
selected cata l ogue o f the COSMOS/UltraVISTA field from 
iMuzzin et al.l ( l20I3bh . The catalogue contains PSF-matched 
photometry in 30 photometric bands covering the wave- 
length range 0.15jum - 24yum and includes the avail- 
able GALEX ([Martinet all I2OO5I). CF HT/Subaru dCapak et al.l 
120071). UltraVISTA (iMcCracken et al.ll2012l) . and S-COSMOS 
(ISanders et al.l l2007h datasets. Source s are selected from the 
DRl UltraVISTA K^-band imaging ( iMcCrackenet al.l l20I2h 
which reaches a depth of Kj = 23.4 at 90% completeness. A 
detailed description of the photometric catalogue construction, 
photometric red shift measurements, a nd stellar mass estimates 
is presented in iMuzzi n et al.l (l20I3bl) . In the next section we 
estimate these properties for the galaxies selected in GCLASS 
in a similar way. In Appendix |B] we show a comparison be- 
tween the UltraVISTA field SMF and the SMF measured in 
GCLASS outside of the clusters. In general the agreement is 
good, even though the GCLASS data are much shallower and 
contain fewer filters. At the low-mass end of the SMF there are 
some small differences due to incompleteness of GCLASS. We 
use UltraVISTA to correct this and provide an unbiased measure 
of the Schechter parameters in the field. 

3. Analysis 

3.1. Photometric redshifts 

We estimate photometric redshifts (photo-z's) using the pub- 
licly a vailable code EAZY dBrammer et al.l2008l) . This code was 
tested dHildebrandt et al.ll20I0h and performs very well on sim- 
ulated and real imaging data. Input to the code are fluxes in the 
10 or 1 1 available bands and their errors. 

We checked for possible systematic problems in the photo- 
metric calibration or photo-z code by comparing the estimated 
photometric redshifts with spectroscopic redshifts where the 
samples overlap, see Fig. [1] The performance is then quantified 
by the scatter, bias and outlier fraction of this comparison. First 
we calculate Az = "'''!"' '"'"'"' for each object with a reliable spec- 
troscopic redshift. For historical reasons and to facilitate com- 
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Fig. 1. Spectroscopic versus photometric redshifts for the 10 
GCLASS clusters. Outliers, objects for which Az > 0.15, are 
marked in red. The outlier fraction is less than 5%, the scatter of 
the remaining objects is cr~ = 0.035. 



parisons with other photo-z studies, we define outliers as objects 
for which |Az| > 0. 15. For the remaining measurements we mea- 
sure the mean of |Az| and the scatter around this mean, cr.. We 
find the following typical values: a scatter of cr~ = 0.035, a bias 
of |Az| = 0.005, and fewer than 5% outliers. Specifically, in the 
redshift range of the clusters (0.867 < z < 1 .335), we find a scat- 
ter of cr- - 0.034, a bias of |Az| < 0.005, and about 7% oudiers. 
We find that the scatter is in the range 0.031 < cr- < 0.042 for 
all 10 clusters, and therefore the differences in photo-z quality 
between the clusters is negligible. 

Whereas these values are computed for the entire population 
of galaxies, a subdivision by galaxy type shows that photo-z esti- 
mates for quiescent galaxies are more precise {cr. = 0.030) than 
for star-forming galaxies (cr, = 0.039) because of the stronger 
4000A feature in the broad-band SEDs of quiescent galaxies, 
and the presence of emission lines and a larger diversity of in- 
trinsic SEDs in the star-forming population. We therefore make 
the separation by galaxy type when correcting for cluster mem- 
bership in Sect. 13.41 The scatter in photo-z estimates increases 
for fainter objects, however we take this effect into account when 
we correct for cluster membership. 

3.2. Stellar masses and completeness 

We estimate stell ar masses for all objects using the SED fit- 
ting code FAST (iKrieket al.ll2009l) . The redshifts are fixed at 
the measured spec-z, whenever available. Otherwise we use the 
photo-z from EAZY, and the stellar population libraries from 
iBruzual & CharlotI (l2003b are used to obtain the model SED that 
gives the best fit to the photometric data. We use a parameter- 
ization of the star formation history as SFR oc e^'^'^, where 
the time-scale r is allowed to ra nge b etween 10 Myr and 10 
Gyr. W e also assume a Chabrieg (l2003b IMF, solar metallicity, 
and the ICalzetti et al.l (l2000i) dust law. These settings are iden- 



tical to those used for the measurement of stellar masses in the 
Ultra VISTA sample, in order to provide a fair comparison. 

We estimate the mass completeness limits for each of the 
clusters in the following way. First we measure the depths of the 
Ks detection bands by performing simulations in which we add 
artificial sources to these images for a range of magnitudes. We 
then run SExtractor with the same settings as for the construc- 
tion of the catalogue (Sect. 12.3.1b . The recovered fraction as a 
function of magnitude for the simulated sources provides an es- 
timate for the depth of the detection image. Note that the clusters 
at higher redshift were prioritized to have longer exposure times 
and therefore deeper detection bands, leading to more homoge- 
neous detection limits in terms of absolute magnitude and stel- 
lar mass. Magnitude values corresponding to the 80% recovery 
limit, which are typically ~ 22.5magAB, are given in Table[T] 

We estimate stellar mass limits that correspond to these 80% 
Ks-band completeness limits in two different ways. The first 
method uses the UltraVISTA catalogue, which is about a magni- 
tude deeper than GCLASS in the Kj-band. For each cluster we 
select all galaxies from UltraVISTA with a photometric redshift 
within 0.05 from the cluster redshift. By comparing the total Kj- 
band magnitudes with estimated stellar masses in this redshift 
range, we identify the galaxy with the highest mass around the 
limiting detection magnitude. This is the object with the high- 
est mass-to-light ratio, corresponding to the reddest galaxy in 
UltraVISTA. All galaxies more massive than these mass limits, 
which are listed in Table [T] will be detected with a probability of 
> 80% in GCLASS. 

Secondly, to provide a comparison, we also give the mass 
limit corresponding to a maximally old stellar population with 
no dust (Bruzual & Chariot 2003), at the redshift of the cluster 
with a flux equal to the detection limit. The mass limits resulting 
from this approach are also given in Table [1] and are similar to 
the first estimates to within several hundredths of a dex for most 
of the clusters. 

Note that for cluster SpARCS-0036 we use the J-band as 
the detection band instead of the Kj-band because the Ks-band 
is of much lower quality. Because the seeing in the J-band is 
significantly better, a J-band selection leads to a stellar mass de- 
tection limit that is 0.3 dex lower than could be obtained with 
a Kj-selection. In Table[T]we therefore give the estimates corre- 
sponding to the J-band. 

3.3. Rest-frame colours 

In the following we make a separate comparison between the 
SMF for star-forming and quiescent galaxi es, and correct eac h 
of the galaxy types fo r cl uster membership. I Wuvts et al.l(l2007h . 
IWilliams et ail(l2009h and lPatelet an(l20I2h have shown that the 
U-, V- and J-band rest-frame fluxes of galaxies can be combined 
into a UVJ diagram to distinguish quiescent galaxies from star- 
forming galaxies, even if the latter population is reddened by 
dust extinction. 

After estimating redshifts for all objects in the photometric 
catalogue, we use EAZY to interpolate the input SED to obtain 
the U-V and V-J rest-frame colours for each galaxy. In Fig. |2] 
we plot those colours for all Kj-band selected objects with 
M* > 10'" Mq. The greyscale distribution shows the galaxies 
in GCLASS that are in the redshift range 0.85 < z < L20, 
but are not part of the clusters, while the red points show the 
objects that are separated from the BCG by less than 2 arcmin, 
and have a photometric redshift within 0.1 from the cluster 
redshift. We select as the quiescent population galaxies with 
(U-V)> 1.3a(V-J)< 1.6a(U-V)>0.88(V-J)h-0.6 
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Fig. 2. Rest-frame U-V versus V-J colours for galaxies with stel- 
lar masses exceeding 10'° M© to differentiate between quiescent 
and star-forming galaxies. The arrow indicates the reddening 
vector for dust. Combining both colours facilitates a distinc- 
tion between both galaxy types, even when there is significant 
reddening by dust. In grey is the distribution of galaxies from 
GCLASS that are between 0.85 < z < 1 .20, but outside the clus- 
ters. A relative density scale is provided. The red points show 
photo-z selected cluster members with projected position less 
than 2 arcmin from the cluster centres. 



(e.g. IWhitaker et all 1201 ih . This dividing line is shown in the 
figure. For reference, the dust-reddening vector is also shown, 
indicative of a dust-independent separation of quiescent and 
star-forming galaxies. 

Comparing the cluster and field galaxy populations, we find 
that 68% of the cluster galaxies in this mass range are quiescent, 
whereas only 42% of the field galaxies are quiescent. This shows 
that the cluster population is dominated by quiescent galaxies, 
whereas the field has a more mixed population of galaxies. Note, 
however, that the distribution of colours due to dust-reddening 
within the separate galaxy types is similar for the two environ- 
ments. 

3.4. Cluster member selection 

Due to the scatter in the photometric redshift estimates, selecting 
cluster galaxies based on photometric redshifts will result in the 
sample being contaminated by fore- and background galaxies. In 
this section we combine the photometric Ks-band selected multi- 
colour catalogue and the sub-sample of galaxies with spectro- 
scopic information to select a complete sample of cluster mem- 
bers. We will use the following terminology. By "false positive" 
we refer to an object that is not part of the cluster, yet has a 
photo-z that is consistent with the cluster redshift. A "false neg- 
ative" is an object that belongs to the cluster, but has a photo-z 
inconsistent with cluster membership. A "secure cluster" object 
is correctly classified as being part of the cluster based on the 
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Fig. 3. This is an adaptation of Fig. [T] and shows a composite 
plot of the 10 clusters to measure the fraction of false positives 
and false negatives, after separating quiescent and star-forming 
galaxies. By plotting the difference with respect to the cluster 
redshift, all clusters are effectively plotted on top of each other. 
This plot shows that the Zphot measurements for star-forming 
galaxies have a larger scatter than the measurements for quies- 
cent galaxies. What is not shown in this plot, is how the purity 
fractions change as a function of mass and radial distance. In 
the analysis we also take account of this mass and radial depen- 
dence; see Fig.|4] 



photo-z, while a "secure field" object is correctly identified as 
being outside of the cluster. 

Given the relatively small fields in which we measure the 
cluster SMF, field-to-field variance complicates a full statisti- 
cal interloper subtraction that is based solely on photometric 
data. However, owing to the extended spectroscopic coverage 
of GCLASS, we can estimate the field contamination from these 
data without having to rely on the statistical subtraction of an 
external field. This way we take account of both false positives 
and false negatives in the photometrically selected sample. The 
objects in the spectroscopic sample were prioritized by 3.6//m 
IRAC flux and prox imity to the cluster core, see Sect. 12.21 and 
iMuzzin et al.l ( |2012|) . This selection leads to a targeting com- 
pleteness that is, to first order, a function of radial distance and 
stellar mass only. 

For these targets we measure the differences between photo- 
z's and the redshift for each cluster, and between spec-z's and the 
redshift of the cluster. A composite for all 10 clusters is shown 
in Fig. [3j after separating between quiescent and star-forming 
galaxies. This can be considered as a different representation of 
Fig. [1] where the data for all clusters have been folded on top 
of each other. Galaxies with |Az| < 0.05 are selected as prelimi- 
nary cluster members based on their photometric redshifts. The 
red crosses show false positives, orange crosses indicate false 
negatives. Green (blue) symbols show objects that are identified 
as secure cluster (field) galaxies. Note that, although we could 
have started with any cut on |Az|, the |Az| < 0.05 criterion con- 
veniently yields a number of false positives that approximately 
equals that of false negatives. 

For the objects in the photometric sample that do not have a 
spectroscopic redshift, we use these fractions of false positives 
and false negatives to correct the number counts for cluster mem- 
bership. To make sure that the spec-z subsample is representa- 
tive of the photo-z sample, we have to estimate this coiTection 
as a function of radial distance and stellar mass. This separa- 
tion ensures that we take account of the spectroscopic targeting 
completeness, the mass- and radially-dependent overdensity of 
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Fig. 4. The correction factors for the photo-z selected members that have no spec-z information, estimated from the subsample of 
objects that do have spectroscopic overlap. A separation by radial distance and stellar mass is made, and these factors are multiplied 
to yield the total correction factor for each galaxy. A correction factor > 1 indicates that the number of false negatives exceeds the 
number of false positives in that bin. In the bottom panels the spectroscopic targeting completeness is shown. 



the cluster compared to the field, and the flux dependence of 
the photo-z quality. In Fig.|4]we show the correction factors, as a 
function of radial distance (left panel) and as a function of stellar 
mass (right panel). Error bars give the 68% confidence regions 
estimated from a series of Monte-Carlo simulations in which 
we randomly draw a number for secure cluster members, false 
positives and false negatives from a Poisson distribution in each 
mass-, and radial bin. A correction factor > 1 indicates that the 
number of false negatives exceeds the number of false positives 
in that bin. Corrections are roughly constant with M*, decreas- 
ing slightly at large radii, but the selection of photo-z members 
as objects with |Az| < 0.05 ensures that the corrections are small 
in general. If we change the cut to 0.03, 0.07 or 0.10, this leads 
to different membership coiTections. However, after these cor- 
rections have been applied, we find that these cuts give results 
that are consistent within the errors. 



Down to the mass-completeness of the clusters there are 283 
spectroscopically confirmed cluster members. We divide the 255 
photo-z members for which we do not have spectra over a 2- 
dimensional array of 3 radial bins and 8 stellar mass bins, and 
correct them for membership by multiplying with both the ra- 
dial and mass-dependent correction factors (as shown in Fig.|4]i. 
Because the corrections are relatively small, the way we bin the 
data only has a minor eff'ect on the results. The dominant source 
of uncertainty is of statistical nature. 



4. Results 

4. 1 . The cluster stellar mass function 

We measure the cluster galaxy SMF from the sample of galaxies 
in the 10 GCLASS clusters, obtained as described in Sect. 13.41 
This is done by summing over the 3 radial bins so that we mea- 
sure the SMF out to a projected radius of 1 Mpc. The summation 
is done separately for quiescent and star- forming galaxies, which 
were identified using the UVJ criterion (Sect. 13.3) . The errors 
from the Monte-Carlo simulations that we discussed in Sect. 13.41 
are propagated. Note however that the spectroscopic targets only 
contribute a Poissonian error, since these do not need to be sta- 
tistically corrected for cluster membership. 

The blue points in the left panel of Fig. |5]show the SMF for 
the star- forming galaxies in the 10 clusters, while the red points 
show the quiescent population in the clusters. The total galaxy 
SMF is the sum of the two galaxy types, and is shown in black. 
The fraction of quiescent and star-forming galaxies to the to- 
tal number of galaxies is shown in the bottom panel. The data 
points are also given in Table [3] Note that the quiescent pop- 
ulation dominates the SMF of the cluster galaxies over almost 
the entire mass range we study. The BCGs are not included in 
this plot, nor in the rest of the analysis in this paper. Although 
the satellites in the galaxy clusters are believed to originate from 
an infalling population of centrals in the field, the BCGs are the 
central galaxies in massive cluster haloes and do not have a field 
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Fig. 5. Comparing the cluster SMF (left panel) with a similar representation of the field SMF (right panel). The total SMFs (black 
points) are separated by galaxy type. Red points show the quiescent galaxies and the blue points show star-forming galaxies. The 
best-fitting Schechter functions are overplotted for each SMF sample. Note that the red points have been offset by 0.01 dex for better 
visibility. In the bottom panel we show the fractional contribution of quiescent and star-forming galaxies to the total population, and 
the curves show the fractional contributions of the Schechter functions. The relative contribution of quiescent galaxies is shown to 
be higher in the cluster than in the field. Note that the error bars on the field data are smaller than the data point symbols, because 
only Poissonian errors are taken into account. 



counterpart. Consequently, BCGs do not follow the Schechter 
function that describes the rest of the cluster galaxi es. For a study 
of the stellar mass evolution of BCGs we refer to iLidman et al.l 

Because the selection bands of some of the clusters are 
not sufficiently deep to probe the SMF down to 10"^ M© (see 
Table[T]i, the lowest two stellar-mass bins are composed of galax- 
ies selected from 6 and 7 clusters, respectively. These two bins 
were scaled up by assuming the richnesses of the clusters are 
similar, i.e. multiplying their values with a factor of ^ and y, 
respectively. A rough estimate of the richnesses of the individual 
galaxy clusters shows that these corrections factors are accurate 
to within 10%. 

We perform a small additional completeness correction 
based on a comparison of the field SMF measured from 
UltraVISTA and the parts of GCLASS that are outside the 
clusters (i.e. the field SMF from GCLASS; see Appendix |B]i. 
Because of the depth of its photometry, UltraVISTA is complete 
at M* > 10"' Mq in the redshift range 0.85 < z < 1.20. We 
compare the field estimates in Appendix |B] and find that there 
is a good quantitative agreement in both the shape and normali- 



sation of the field SMF between the surveys at this stellar mass 
range, except for the lowest three mass points. This suggests that 
there may be residual incompleteness in GCLASS that affects 
the lowest mass points. Assuming that this incompleteness af- 
fects the cluster and field data of GCLASS in a similar way, we 
correct the GCLASS cluster SMF points for the star-forming and 
quiescent galaxies with small factors, up to 37% at the lowest 
mass bin for the quiescent galaxies. This correction changes the 
best-fit Schechter parameters for the cluster fits in the following 
way. M* increases by 0.01, 0.10 and 0.08 dex and a becomes 
more negative by 0.07, 0.33 and 0.26 for the total, star-forming 
and quiescent population, respectively. These changes do not af- 
fect any of the qualitative results in this paper, nor change the 

conclusion in any way. 

We fit a Schechter dSchechterll 19761) function to the binned 
data points. This function is parameterized as 



(l)(M) = ln(10)O*[l0 



(M-M')(l+a) 



]-exp[-10<*^-^'']. 



(3) 



with M* being the characteristic mass, a the low-mass slope, and 
d>* the total overall normalisation. 
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Table 3. The values for the data points of the galaxy SMF that are shown in Fig. [5] These are the raw, membership-corrected, 
numbers of galaxies for the clusters. To obtain the units shown in the figures for the clusters, these values need to be multiplied by 
5, since the binsize is 0.2 dex in stellar mass. Numbers in brackets show the total number of spectroscopic cluster members in each 
bin. Note that the spectroscopic completeness is highest in the high-mass bins. Errors represent Icr uncertainties estimated from 
Monte-Carlo simulations for the cluster data, and Poissonian errors for the field data. 
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Because the number of sources in the brighter stellar mass 
bins is low, we are in the regime where errors cannot be rep- 
resented by a Gaussian distribution and therefore ordinary x^ 
minimisation is not appropriate. Consequently, we take a gen- 
eral maximum likelihood approach where we use the probability 
functions on each data point obtained from the Monte-Carlo sim- 
ulations. This way we compute the likelihood function X on a 3 
dimensional grid of Schechter parameters. The best fitting pa- 
rameters M* and a, corresponding to Xmax, are listed in Table |4] 
and the corresponding Schechter function is shown in the left 
panel of Fig. |5] (black curve). The Schechter function provides 
a reasonable fit to the data, with a Goodness of Fit (GoF) of 
1.82. We also give the 68.3% confidence levels on each param- 
eter after marginalising over the other two parameters. We take 
this confidence interval to be the region where 2 ln(Xmax/-C) < 1 . 
However, since these parameters are known to be degenerate, we 
show confidence regions in Fig. Q The black curves correspond 
to 68.3% and 95.4% confidence levels after marginalising only 
over 0*. 

In general, uncertainties on individual mass measurements 
of the galaxies lead to a bias in th e shape of the SMF and the 
best f itting Schechter parameters (lEddingtonlll913l: iTeerikorpil 
l2004l) . Especially for high masses, where the slope of the SMF 
is steep, the shape of the SMF can be biased because of galaxies 
scattering to adjacent bins. To study the magnitude of this effect 
on our analysis, we need to estimate the stellar mass scatter of 
galaxies in each bin of the SMF. To do this we created 100 Monte 
Carlo realisations of the photometric catalogue, in which we ran- 
domly perturb the aperture fluxes following the estimated statis- 
tical errors on these measurements. Then we estimate photo-z's 
and stellar masses for the entries of these catalogues in a similar 
way as for the standard analysis. At the high-mass end, where 
the SED fitting is mostly supported by spec-z's (see Fig. |4]or 
Table |3]l, the scatter is about 0.05 dex in stellar mass. For lower 
masses the scatter increases towards 0.08 (0.10) dex in stellar 
mass for quiescent (star-forming) galaxies. Even if all galaxies 
would scatter to higher masses, the bias in the Schechter parame- 
ter M* would be 0.05 dex. I n reaUty a might also cha nge slightly 
due to Eddington bias (e.g. Ivan der Burg et al.ll2010l) . but we ex- 
pect the bias of the combination of Schechter parameters to be 
substantially smaller than the size of the l-cr statistical error con- 
tours in Fig. [T] Given also that the systematic uncertainties due 
to assumptions regarding the IMF, models on the stellar popu- 
lations, star-formation histories and metallicity, are substantially 



Table 4. A comparison between the best fitting Schechter param- 
eters and their 68% confidence intervals for the different galaxy 
types and environments. 
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Goodness of Fit (GoF) defined as ^-/dof for the field survey (we 
assumed Gaussian statistics owing to the large number counts in 
this survey). For the cluster fits we used an analogous expression 
from the Maximum likelihood fitting method. 



larger (e.g. lMarchesini et al.ll2009|) . we do not attempt to correct 
for this bias in the current analysis. 

4.2. Cluster versus Field 

We compare the cluster results with the field galaxy SMF by 
selecting all galaxies with a photometric redshift in the range 
0.85 < z < 1.20 from the Ultra VISTA survey. Since the 
Ultra VISTA survey is superior in depth compared to GCLASS, 
the SMF can be measured down to 10"^ M© in this redshift range. 
The right panel of Fig. |5] shows the field total SMF in black, 
which is composed of 13633 galaxies in this mass and redshift 
range. The best fitting Schechter function for the field sample is 
found by minimizing ;^'- on a 3 dimensional grid of Schechter 
parameters, and is represented by the black curve in the right 
panel of Fig. |5] For a comprehensive comparison between the 
S MF from Ultra VIS TA and other field estimates we refer to 
Muzzin et al.l ( l20I3ah . There it is shown that the SMF of the en- 
tire galaxy population, measured with this catalogue, is in good 
agreement with previous measurements. 

To better compare the shape of the total SMF in the two en- 
vironments, we refer to the left panel of Fig. |6] where the ma- 
genta points show the galaxy SMF from UltraVISTA, and the 
black points show the SMF for the cluster galaxies. The field 
data have been scaled such that the Schechter functions of the 
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Fig. 6. Galaxy SMFs for different galaxy types and environments. The left panel shows the total galaxy population in the cluster 
(black) and the field (magenta). The middle panel shows the cluster and field SMF for the subset of star-forming galaxies. The 
right panel shows the subset of quiescent galaxies. The field data have been scaled vertically to match the cluster SMF at M* of the 
cluster Error bars show the 68% confidence regions from Monte-Carlo simulations (on the cluster data), or Poisson error bars (field 
data). 
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Fig. 7. The 68% and 95% likelihood contours for the Schechter parameters M* and a, after marginalising over the (p* parameter 
Black lines show the cluster contours, while magenta lines show the contours for the field data. -H-signs show the single best fit 
Schechter parameters. The regions corresponding to the cluster SMF were obtained using maximum-likelihood fitting of the Monte- 
Carlo simulated data. 



cluster and field intersect at the characteristic mass M* of the 
cluster The best fitting values for the a and M* parameters are 
given in Table H) with their 68.3% confidence levels. Because 
we only included Poissonian errors on the field SMF data, the 
GoF of the Schechter fits are rather high (up to 4.66 for the total 
galaxy population). At this level of detail it is also possible that 
the Schechter function is no longer an adequate description of 
the data. The magenta contours in the left panel of Fig. |7] show 
the 2-d confidence contours for the field. 



4.3. Star-forming vs Quiescent galaxies 

We separate the Ultra VISTA galaxy catalogue between quies- 
cent galaxies and star-forming galaxies by using their estimated 
rest-frame U- V and V-J colours, as was analogously done for the 



cluster galaxies in Sect. 13.31 We compare the shapes of the SMF 
for each galaxy type between the different environments. 

In the middle panel of Fig.|6]we show the shape of the SMF 
for star-forming galaxies in the field (magenta) and the cluster 
(black), together with their best-fitting Schechter functions. The 
field data have been normalised so that the Schechter functions 
intersect at the characteristic mass M* for star-forming galaxies 
in the cluster The corresponding 68% and 95% confidence re- 
gions for the Schechter parameters a and M* are shown in the 
middle panel of Fig. [T] The best-fitting Schechter parameters and 
their error bars are also given in Table 

4.4. Normalisation of the SMF 

The data points in Fig. |6] are arbitrarily normalised to provide 
for an easy comparison of the shapes of the SMF between the 
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field and cluster samples. As a consequence, the 0* parameters 
corresponding to the best fitting Schechter function have no di- 
rect meaning. Normalised by volume the cluster is, by defini- 
tion, substantially overdense compared to the field. To be able 
to better interpret the differences of the SMF between the field 
and cluster environment in Sect. |5] we therefore normalise the 
SMF by the total amount of matter in each respective part of the 
Universe. 

For the UltraVISTA field reference we take the total comov- 
ing volume within a redshift range 0.85 < z < 1 .20 and an un - 
masked survey area of 1 .62 square degree (iMuzzin et al.l2013bh . 
After multiplying the volume corresponding to this area in this 
redshift range, 5.9 ■ lO^Mpc^, by the average matter density of 



All galaxies 



the Universe, being 2.8 -10 g cm"^ in our cosmology, we find 
that the total amount (i.e. dark matter -i- baryonic) of matter in 
this volume is about 2.4 ■ 10'^ Mq. 

Given the values for M200, which are presented in Table |2] 
and are based on the dynamical analysis of the GCLASS spec- 
tra, we estimate t he concentration par ameter corresponding to 
the NFW p r ofiles (iNavarro et al.lll996l) for these systems from 
iDuffv et all (l2008h . We integrate these NFW profiles along the 
LOS and out to a projected radius of 1 Mpc , yielding a total mass 



of 5.6 X 10'^ Mq for the 1 clusters. Since ISheldon et al.l (l2009l) 



and lHoekstra et al.l (l2000h have shown that, although cluster cen- 
tres are dominated by luminous matter, the mass to light ratio 
(M/L) of clusters within a distance of 1 Mpc is similar to the 
cluster M/L within larger distances, this ensures that we mea- 
sure and normalise the SMF in a representative volume. 

Fig.[8]shows the total SMF for the cluster and the field, after 
normalising by the total masses estimated above. Note that there 
is, per unit total mass, a strong overdensity of galaxies at all 
stellar masses we probe in the cluster environments. In the stellar 
mass range we study, the overdensity ranges from a minimum 
factor of 1.4 at 10'" Mq to a maximum factor of 3.1 at 10"' Mq. 
This shows that the clusters contain a very biased population of 
galaxies, where a relatively high fraction of the total baryonic 
mass is transformed into stars. The field, in contrast, contains 
regions such as voids, where the star formation efficiency is very 
low. 



5. Discussion 

In this section we discuss the implications of the results from 
Sect. |4] We discuss in Sect. ED the shape of the SMF for star- 
forming galaxies, quiescent galaxies, and the total galaxy popu- 
lation. We make a comparison between the cluster and field, and 
also compare our results to measurem ents from the litera ture. We 
proceed to apply a simple model that lPeng et al.l (l2OI0l) showed 
to give a good fit t o the SMF measured at z = from SDSS data. 
iPeng et aD (l2010h could not explore the area of high-z clusters 
with COSMOS and SDSS data, so we confront our results at 
z = 1 with the predictions of their model. 



5. 1 . The shape of the galaxy SMF 

5.1.1. Star-forming galaxies 

Fig. |6] shows that the shape of the galaxy SMF for the sub- 
set of UVJ-selected star-forming galaxies is similar between 
the clusters from GCLASS and the field from UltraVISTA. 
Quantitatively, Fig. [7] indicates that the combination of best- 
fitting Schechter parameters differs by about Icr. The low-mass 
slope a is -1.35;^Q2g for the cluster versus -I.IB^qq^ in the field. 
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Fig. 8. Same as the left panel of Fig. |6] but normalised by the 
total mass (dark matter -1- baryonic) in the field sample (magenta) 
and cluster sample (black). Per unit of total mass the cluster has a 
clear overdensity at all stellar masses we probe. Error bars show 
the 68% confidence regions from Monte-Carlo simulations (on 
the cluster data), or Poisson error bars (field data). In the left 
panel of Fig. |6] we provided an easier comparison of the shapes 
of the two SMFs. 



The characteristic mass M* is 10.86+^-^^ ^^'^ 10-65!ooi ^°^ *^ 
cluster and field, respectively. 

We do not make a quantitative comparison between the liter- 
ature and our measurements of the SMF for star-forming galax- 
ies because the way these star-forming samples are selected is 
different for different studies. Whereas we select a subset of 
star-forming galaxies based on the UVJ-diagram, most other 
studies use either a single colour or a morphological selection. 
Nonetheless, the finding that the shape of the star-forming SMF 
is independent of environment is qualitativel y consistent with 
lower redshift measurements presented by e.g. iBolzonella et al.l 
(l2010h . Note however that the clusters in GCLASS constitute 
much higher overdensit ies than the hig hest densities in the 
COSMOS fields used bv IBolzonella et all (|2010|) . The shape of 
the star-forming galaxy SMF is also measured to be roughly 
constant with cosmic time (e.g. lllbert et al.l2()10HBrammer et all 
l2011h . This shows that, whatever processes are responsible for 
the quenching of star formation in galaxies, they have to operate 
in such a way that the SMF of star-forming galaxies does not 
change shape, even in the highest density env i ronme nts. This is 
a fundamental assumption for the iPeng et all (1201 Ol) quenching 
model that we employ in Sect. 



5.1.2. Quiescent galaxies 

Fig.|6]shows that for the selection of quiescent galaxies based on 
the UVJ criterion, the shape of the SMF for those galaxies is also 
similar in the different environments probed by GCLASS and 
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Ultra VISTA. The best fitting a for the clusters is -0.28+°-^^ ver- 
sus -0.43;^Q Q4 in the field. Given the degeneracy between a and 
M* , the combination of these Schechter parameters, as shown 
in Fig. |7] also agrees to better than Icr between the field and 
cluster. It seems remarkable that, whatever quenching processes 
are responsible for the build-up of the quiescent population in 
these contrasting environments, they work in such a way that 
the resulting SMF for quiescent galaxies at M* > 10'" M© has a 
similar shape in both environments. 

iRudnick et al.l (|2009|) measured the cluster galaxy luminos- 
ity function of red sequence galaxies in the redshift range 0.4 < 
z < 0.8 and compared their measurements with the field lumi- 
nosity function. They also found little difference in the shape 
of the quiescent lumi n osity function between the two environ- 
ments. iRudnick et alJ (!2009) also found a hint of a shallower 
low-mass slope in the cluster compared to the field. Note that 
they use a different selection of red galaxies, so that their red 
sequence selected sample might be contaminated by reddened 
star-forming galaxies. 



5.1 .3. The total galaxy population 

Whereas the SMF for each of the galaxy types appears to be 
similar in the different environments probed by GCLASS and 
UltraVISTA, Figs.|6]&|2]show that the SMF for the total galaxy 
population is significantly different. This is because the frac- 
tion of quiescent galaxies is higher in the cluster. That makes 
the low-mass slope of the total SMF shallower in the cluster 
compared to the field (see Fig. |6|l. This result is also consistent 
with the measurements s hown for more moderate overdensities 
in the COSMOS field bv iBolzonella et all (l2010l) . We compare 
our results to the literature results from the WIN GS, ICBS and 
EDisCS clusters probed in IVulcani et al.l (l2013h . although our 
sample is unique in this combination of redshift range and pho- 
tometric depth. 

IVulcanietaT1(l2013h assumed a Kroupa (lKroupal200lh IMF, 
which yields stellar masses consistent with Chabrier to within 
several 0.01 dex. For the sample of WINGS clusters (0.04 < 
z < 0.07) they measure Schechter parameters M* = 10.82 + 
0.13 and a - -0.88 + 0.31. Although the redshift distribution 
is very different from the GCLASS sample, they agree within 
1 -cr with the contours shown in Fig.|7] The measured Schechter 
parameters for the ICBS clusters (0.3 < z < 0.45) are M* - 
11.37 + 0.28 and a = -1.29 + 0.41. Note that this point lies in 
the direction of the correlation between M* and a, as is shown 
in Fig|7] The same is true for the EDisCS clusters (0.4 < z < 
0.8), for which rVulcani et all (1201 3h report Schechter parameters 
M* = 11.15 + 0.07 and 0- = -1.03 + 0.08. 

Another fundamental observable of a population of galaxies, 
besides their SMF, is the d i stribu tion of specific star formation 
rates (sSFRs). IWetzel et al.l(l2012h studied the distribution of sS- 
FRs for central and satellite galaxies as a function of stellar mass 
in a range of environments. They show that the distribution of 
sSFRs is clearly bimodal, with clear populations of active and 
passive galaxies. Interestingly, they show that the location and 
shape for each of the two peaks is independent of environment, 
and that only the relative amounts of star-forming and quiescent 
galaxies o ccupying the peaks d iffer as a function of environment. 
Likewise. iMuzzinet all (|2012|) show that for the GCLASS data 
the sSFR of star-forming galaxies in a given mass bin is also 
independent of environment. These results are analogous to our 
measurements for the SMF, which can also be considered as a 
sum of the quiescent galaxy SMF and the star-forming galaxy 



SMF. Having a different fraction of quiescent galaxies in op- 
posing environments, the total galaxy SMF will look different 
wher eas the SMF for ea ch galaxy type is similar, analogous to 
what I Wetzel et all (1201 2h found for the distribution of sSFRs. 

5.2. A simple quenciiing model 

It has been known for several decades that the fraction of quies- 
cent galaxies increases with both stellar mass and environmental 
density (e. g. Baldrv et aL 2006h. How ever, recent studies (e.g. 
Pe ng et al.l 120 10: Muz zinet alTl2012h have suggested that the 
quenching of star-forming galaxies can be fully separated in two 
distinct quenching tracks, dubbed "mass quenching" and "en- 
vironmental quenching". The assumption that the shape of the 
galaxy SMF for star-forming galaxies is universal, which is sup- 
ported by our measurements, places constraints on the way these 
quenching processes operate. 

To interpret our observed d ata in this c ontext we consider 
the simple model proposed by iPeng et all ( |2010|) . This model 
is based on the observed constancy i n the shape of the SMF 
for star-forming galaxies with redshift. iPeng et aD (12010) use a 
combination of mass quenching and environmental quenching, 
processes which they assume to act independently of each other, 
to build up the quiescent population. The basic descriptions for 
these quenching tracks are demanded to operate such that the 
shape of the SMF for star-forming galaxies is independent of 
environment. 

Because star-forming galaxies are forming stars at a rate 
that scales roughly linearly with their stellar mass (the ob- 
served sSFR for this population is roughly independent of mass 
dNoeske et al.ll2007h X mass quenching is supposed to preferen- 
tially quench high mass galaxies in order to keep the SMF for 
star-forming galaxies unchanged. Therefore the resulting galaxy 
SMF for this quenched population is expected to contain an 
excess of high mass galaxies and hence has a shallow low- 
mass slope. In high-density environments the fraction of quies- 
cent galaxies incr eases compared to low-density environments. 
iPeng et all ( 1201 Ol) assume that this increase is caused by the pro- 
cess of environmental quenching. The environmental quenching 
efficiency is assumed to be independent of mass, so that the re- 
sulting SMF of the environmentally-quenched galaxies has the 
same shape as the star-forming galaxy SMF. With some addi- 
tional quenching d ue to, what they presume to be, merger pro- 
cesses, iPeng et al.i ( 1201 Ol) showed that this model works very 
well at reproducing the SMF mea sured in the redshift range 
0.02 <z< 0.085 from SDSS DR7 dAbazaiian et al.ll2009l) data. 
The regime of z ~ 1 clusters however has not yet been tested 
against their model, since this range is not probed by COSMOS . 

The model, however, makes predictions for the SMF 
at higher redshif ts over a range of environmental densities 
( iPeng et al1l2010l Fig. 13), and we compare these predictions 
at z = 1 to the SMFs measured from GCLASS and UltraVISTA. 
The predictions from their model are separated by environmen- 
tal density in four quartiles, with Dl (D4) corresponding to the 
lowest (highest) density quartile. The (especially environmen- 
tally) quenched galaxies contribute more substantially to the to- 
tal galaxy population in D4 compared to Dl, leading to a higher 
fraction of quiescent galaxies. The left panel of Fig.|9]compares 
the prediction of the highest environmental density quartile (D4) 
with the measurement of the cluster galaxy SMF from GCLASS. 
We fitted the total normalisation as a free parameter, but left the 
relative normalisations of star- forming and quiescent galaxies 
unchanged. Note that the GCLASS clusters constitute the most 
massive structures around z = 1 and therefore contain higher 
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Fig. 9. The left panel compares the lPeng et al.l (1201 Ol) model prediction in the environmental density quartile D4 with our GCLASS 
cluster SMF measurements, which were already presented in Fig.|5]and Table[3] A separation between the two quenching processes 
is made. The right panel makes a similar comparison between the model in quartile Dl with the UltraVISTA field data. Note that 
the relative normalisations of star-forming and quiescent populations are fixed, and that the populations are fitted simultaneously. 



overdensities than the D4 reference. Nevertheless, the D4 model 
provides a reasonable fit to the data, where the quiescent fraction 
of galaxies between the model and the data is well matched. In 
future studies it would be interesting to compare the iPeng et all 
(1201 ol) model for the upper 5% in environmental density to the 
cluster data, which would be a closer match to their density. 

The UltraVISTA field is expected to contain a range of envi- 
ronmental densities. The measured SMF from these data should 
therefore be matched to a combination of the D1-D4 models. 
However, the right panel of Fig. |9] shows that even the low- 
est environmental density quartile Dl overpredicts the quies- 
cent fraction of galaxies in the field of UltraVISTA at z = 1 . 
The caveat is that the separation of star-forming and quiescent 
galaxies is done differently between the data and the model. 
Our sample of UVJ-selected star-forming galaxies includes star- 
forming galaxies that are seen edge-on and therefore reddened 
b y dust, whereas a r est-frame (U-B) colour selection, as applied 
in lPeng et al.l (1201 Ol) . identifies these objects as being on the red 
sequence. 

To reconcil e the ap parent disagreement between the data and 
the iPeng et all (1201 Ol) model in predicting the quiescent frac- 
tion of galaxies, we consider the following simplified analytical 
model where we only assume mass quenching and environmen- 
tal quenching, and no additional merging. W e apply this simpli - 
fied model, based on the same principles as i Peng et alj (120101) . 
to the GCLASS cluster data, but use the Schechter fits to the 
UltraVISTA data as a starting point. UltraVISTA is the limiting 
case where the dominant quenching process is mass-quenching. 

We fit the cluster data by a combination of three func- 
tions that describe populations of star-forming galaxies, mass- 
quenched quiescent galaxies and environmentally-quenched qui- 
escent galaxies. Two of these functions are given by the 
Schechter fits that were measured for the UltraVISTA field pop- 
ulation. The quiescent population of UltraVISTA is expected to 
be primarily mass-quenched at the stellar mass range we study. 



so we take the shape of this SMF for the mass-quenched popu- 
lation and allow the normalisation to be adjusted by the fit. The 
SMF for star-forming galaxies is also taken from UltraVISTA, 
and since the functional form of this distribution is assumed to 
be independent of environment, we use the shape of this SMF 
and allow for a change in the normalisation. The third SMF, that 
describes the population of environmentally-quenched galaxies, 
is assumed to have the same shape as the SMF of star-forming 
galaxies, but the normalisation can be adapted in the fit. The 
sum of the functions for both quenched populations is fitted to 
the data points that describe the SMF for quiescent galaxies. 

Now that the functional forms of the three populations that 
we fit are defined, the normalisations are adapted by fitting two 
free parameters in the following way. One free parameter x de- 
scribes how the three functions move relative to each other, and 
constrains the percentage of star-forming galaxies that is envi- 
ronmentally quenched by the cluster. This gives rise to a popula- 
tion of environmentally-quenched galaxies with a normalisation 
of X compared to the star-forming galaxies. The star-forming 
galaxies are reduced by a factor of (1 - x). We do not change 
the relative amount of mass-quenched galaxies. The second free 
parameter describes the total normalisation of these three func- 
tions and has no direct meaning because the cluster and field are 
arbitrarily normalised with respect to each other. 

We perform a maximum-likelihood fit to the data points for 
the star-forming and quiescent galaxies simultaneously, where 
we adapt these two parameters, and find a best fitting value of 
X - 0.44^yy^. Assuming this simple picture we therefore find 
that, besides the quenching processes that also happen in the 
field, the cluster environment has to quench an additional 44% 
of the galaxies to yield the best fit. In the Peng D4 model this en- 
vironmental quenching fraction ranges from 0.17 at 10'" M© to 
0.22 at 10" Mq. Fig.[TO]shows the best fit to the observed SMF 
in the cluster based on this simple model. The blue and red solid 
lines give the simultaneous best fit to both galaxy types, for the 
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star-forming and quiescent populations respectively. The red line 
is composed of a mass-quenched population (dotted red line), 
and an environmentally-quenched population (dashed line). The 
quiescent population at high (> 10'"~ Mq) masses is dominated 
by mass-quenched galaxies, while the population at lower stellar 
masses is dominated by environmentally-quenched galaxies. 
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Fig. 10. The cluster data, which were already given in Table |3] 
with the predictio ns from a simple quenching model based on 
iPeng et al] (1201 Ol) in its most basic form, using Ultra VISTA as 
a starting point. Mass quenching and environmental quenching 
are assumed to act independently. The blue Schechter function 
is the best fit to the star-forming galaxy population and the red 
solid line gives the best fit to the quiescent population. The 
red Une is composed of a mass quenched population (dotted 
red line), and an environmentally quenched population (dashed 
line). This model needs 44% additional environmental quench- 
ing compared to the field to yield the best fit to the data. Note that 
the red points have been offset by 0.01 dex for better visibiUty. 



The best-fitting model does not yield a perfect represen- 
tation of the data, since the model significantly overpredicts 
the number of quiescent galaxies in the low mass regime (< 
10'"'^Mo). At intermediate masses around 10"^-^ Mq the model 
predicts about 30% less galaxies than the data show. The over- 
all Goodness-of-F it for this model is 2. 1 per degree of freedom. 
iPeng et aD (1201 Ol) acknowledge that another term, due to merger 
quenching, is required to fit the data in SDSS and zCOSMOS. 

We know that m ergers occur in clusters (e.g. 
Ivan Dokkum et al] Il999l) . and that the intra-cluster light 
builds up over time, probably by disruptions of relatively low 
mass galaxies (iMartel et al.ll2012h . Also we know that BCG s 
have to grow in stellar mass over time (e.g. iLidman et al.ll2012l) . 
likely by consuming infalling galaxies. It is possible that these 
merging processes are required to reconcile the disagreement 
between the data and this model. 



The intriguing similarity in the shape of the quiescent SMF 
between the cluster and field environments at z ~ 1 su ggests that 
there might be a simpler explanation than the iPeng et al. (201^) 
model that does not involve a large amount of mergers. A simi- 
lar internally-driven quenching mechanism might be responsible 
for the build-up quiescent population in both environments. We 
know that the age of a quiescent galaxy at a given ste llar mass 
does not significantly dep end on its environment (iThomas et al.l 
l2010l:lMuzzin et al]|2012l) . However, for galaxies at a given stel- 
lar mass, their underlying dark matter (sub-)haloes at the time of 
formation might be diff'erent between the clusters and the field. 
"Environment quenching" could therefore refer to an internally 
driven process that is accelerated in cluster sub-haloes compared 
to the field. The finding that the cluster environment has already 
formed a large stellar mass content hy z ~ I (see Fig. [SJ com- 
pared to the field, and that the fraction of quiescent galaxies is 
higher in the cluster than in the field, could be caused by a dif- 
ferent evolution of the underlying dark matter haloes. 

A detailed study of the evolution of the (sub-)halo mass func- 
tion, compared between cluster and field, is required to look into 
the different quenching scenarios. It is required to trace back the 
(sub-)haloes that host the galaxies we study to investigate how 
their progenitors merge and accrete during their formation his- 
tory. Such a study could be useful to better understand the pro- 
cess of environmental quenching. 

6. Summary and Conclusions 

In this paper we measured and compared the galaxy SMF at 
z ~ 1 in the high-density environments probed by GCLASS and 
the field environment from Ultra VISTA. The GCLASS sample is 
composed of 10 rich, red-sequence selected clusters in the red- 
shift range 0.86 < z < 1.34. The Ks-band selected catalogue 
based on observations in 10 - 11 photometric filters allowed us 
to estimate photometric redshifts and stellar masses for galaxies 
in the studied redshift range down to a stellar mass of lO'^M©. 
The extensive spectroscopic sample of GCLASS covers the ma- 
jority of the cluster members, and is critical to account for con- 
taminants in the sample for which we only have photometric red- 
shifts. Galaxies were separated by type (star-forming versus qui- 
escent) based on their rest-frame U-V and V-J colours. For each 
galaxy in the photometric sample we estimated the probability 
that it is part of the cluster based on its type, stellar mass and 
clustercentric distance. This resulted in a statistically complete 
sample of cluster members to measure the SMF from. 

As a reference field SMF we used UltraVISTA, which 
is a new NIR survey that overlaps with COSMOS, resulting 
in 30 band photometric coverage over 1.62 square degrees. 
Analogously to GCLASS, sources were selected from the Kj- 
band, and galaxies were separated between the star-forming and 
quiescent type using the rest-frame UVJ fluxes. This led to a 
measurement of the SMF for field galaxies at 0.85 < z < 1.20 
that is complete down to stellar masses of IO'^Mq. 

We found that the shape of the SMF for star-forming galax- 
ies is similar between the cluster and field environment, and that 
the combination of best-fitting Schechter parameters a and M* 
agree to Icr between the cluster and field. Furthermore, for the 
samples of quiescent galaxies we obtain a similar result. The 
shape of the SMF for quiescent galaxies is similar between the 
cluster and field at M* > 10'°Mq. The shape of the SMF for 
the total galaxy population is significantly different between the 
cluster and field. This is caused by a different fraction of quies- 
cent galaxies in both environments. We find that there is a rela- 
tive deficit of galaxies with low stellar masses in the cluster com- 
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pared to the field. However, when we normalise the SMF by the 
total amount of matter in each respective part of the Universe, we 
find that there is a strong excess of galaxies over the entire stellar 
mass range we probe. This indicates that the cluster environment 
must have been substantially more efficient in transforming mass 
into stars compared to the field. Note that this does not imply that 
field galaxies are less efficient, but rather it is the consequence of 
the fact that voids contain dark matter, but relatively few stars. 

The similarity in the shape of the quiescent and star-forming 
SMF between the cluster and the field indicates that, if different 
processes are to be responsible for the quenching of star for- 
mation in different environments, these processes have to work 
in such a way that the shapes of the quiescent and star-forming 
SMF are similar in these different environments at z = 1 . This 
poses a challenge to analytical models that attempt to explain 
the build-up of the quiescent population by a combination of 
mass quenching and environment quenching. A simple model 
suggests that 44^^% of the star- forming galaxies which normally 
would be forming stars in the field, would be quenched by the 
cluster Although the physical processes that cause galaxies to 
quench environmentally are not yet completely understood, it is 
clear that a process like environmental quenching at z ~ 1 is 
important. 
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Appendix A: Data processing and catalogue 
creation 

This Appendix is meant to give a more elaborate description of 
the data reduction steps (Sect. lA.H and in particular the pro- 
cedure for homogenising the PSF and measuring colours using 
Gaussian weighted apertures (Sect. IA.2b . Because we combine 
photometric data over a wide range of wavelengths and for clus- 
ters that are both in the Northern and Southern sky, we necessar- 
ily have to combine data from different telescopes and/or instru- 
ments. 



A.1. Photometric data reduction 

The standard reduction steps include bias and flatfield correc- 
tions. Although the images have been flatfielded (e.g. by Elixir 
for the MegaCam data) to yield a uniform zeropoint for the 
source fluxes, there are still residual background patterns due to 
scattered light, fringe residuals, and amplifier drift (Cuillandre, 
private communication). These patterns are reasonably stable 
over time, and since most exposures in a given filter have been 
taken consecutively on the same night, we can subtract these 
background effects. We do this by using the dithered pattern of 
observations to differentiate signals that are on a fixed position 
on the ccd array from sky-bound signals. 

To remove cosmic rays from ccd images one usually com- 
pares different frames of the same part of the sky. However, since 
we only have a few deep exposures in some of the filters, the 
number of overlapping frames of our data set is not always suffi- 
cient to be able to identify all cosmic rays. Therefore we remove 
cosmic rays by using the L aplacian Cosmic Ray Identification 
method (Ivan Dokkurrj2001l) . which works on individual images. 
We optimise the parameters in the setup of the code such that we 
minimise the amount of false positives (bright stars) and false 
negatives. We do this by testing the code on a range of images 
with different seeing. The only parameter that has a significant 
influence on the fraction of false positives and false negatives is 
objlim, which we take to be 3.0. 

We obtain astrometri c and relative photometric solutions for 
each chip using SCAMP jBertinl 120061) . where we use all expo- 
sures in a given filter for all clusters at once to effectively in- 
crease the source density, and obtain stable solutions. As a ref- 
erence catalogue we use SDSS-DR7 data, or the USNO-B cata- 
logue whenever a cluster falls outside the SDSS footprint. This 
leads to consistent astrometric solutions between the different 
filters with an internal scatter of about 0.05". 



A.2. PSF liomogenisation and colour measurements 

Because the shape and size of the point spread function (PSF) 
are different between exposures and filters, it is non-trivial to 
measure accurate colours of a galaxy. The simplest approach 
would be to take the ratio of the total flux of a galaxy in dif- 
ferent bands, but this requires very large photometric apertures: 
for background-limited observations these are very noisy. 

However, a reliable colour measurement for the purpose of 
photometric redshift determination can also be made by tak- 
ing the ratio of aperture fluxes in different bands, provided 



these apertures represent the same intrinsic part of the galaxy. 
We have followed this approach here, based on a modification 
of the Gaussia n-aperture-and-PSF (GaaP) photometry method 
(lKuiikeiJl2008h . 

The first step is to convolve each image with a suitable 
position-dependent kernel that modifies its PSF into a uniform 
size, circular and Gaussian. This kernel can be constructed us- 
ing the shapelet jRefregier 2003 ; Kuiiken 2006i) formalism, as 
was d one in the local approach described in iHildebrandt et al.l 
(1201 2h . with one modification: here we allow the resulting PSF 
size for each image to be different. Specifically, for each filter 
and field we choose the size of the resulting PSF to be slightly 
larger (ca. 10%) than the median gaussian radius of all bright 
stars in the images. To obtain a stable PSF in the stacked images 
for each filter we Gaussianise the PSF of the individual astro- 
metrically cor rected exposures before stacking. 

Following KuijkenI (l2008h we then measure fluxes in the 
following way. Instead of using a function where the weight 
is either or 1, as is the case for regular aperture photome- 
try measured with a top-hat weighting function (e.g. by run- 
ning SExtractor in dual image mode), we use a smooth weight 
function that makes use of the fact that the S/N for each pixel 
decreases away from the peak pixel. When the PSF in each fil- 
ter follows a Gaussian profile, the choice to perform photometry 
using a Gaussian weight function is computationally convenient, 
as we show n ext. 

iKuiikenI (l2008h defines the "Gaussian-aperture-and-PSF" 
flux Fq as the Gaussian weighted flux a source would have if 
it were observed with a Gaussian PSF with the same width q as 
the weight function. Hence 



s fdre '■'^2''' rdr'5(r')' 



-(r-r')V2?- 

27Tq^ 



(A.1) 



where S is the intrinsic light distribution of the source (i.e. before 
smearing with the PSF) and q is the scale radius of the weight 
function. It is straightforward to simplify Eq. lA. 1 I to 



H 



5(r)e 



-r-/4g^ 



(A.2) 



which shows that Fq is a Gaussian-aperture photometric mea- 
surement of the intrinsic galaxy. 

After gaussianising the images, S has already been con- 
volved with a Gaussian that has a constant dispersion gpsF for 
each stacked image. The flux distribution on the ccd is therefore 



/(r) 



-(r-r')-/2g>sf. 

2nglsp 



= rdr'5(r')- 
Analytically we have an identical expression for Fg 

J -^^ ~ S PS F 



(A.3) 



(A.4) 



which thus shows that the same intrinsic aperture flux Fg can 
be measured from images with a range of Gaussian PSF sizes. 
Therefore, from our PSF-gaussianised images, we can mea- 
sure colours of the same intrinsic part of the galaxy if we use 
Gaussian weight functions to measure fluxes. Note that it is no 
longer necessary that the stacks of the different filters have a PSF 
with the same Gaussian FWHM, as long as the weight function 
is adapted accordingly for each filter 

We adjust q to ensure the aperture roughly matches each 
galaxy's size, to optimise the S/N. We base our choice for q 
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Table A.l. The GCLASS photometric data set. The instruments used for the different clusters and filters are indicated. The limiting 
magnitudes reported are median 5-cr flux measurement limits for point sources measured with a Gaussian weight function. 



Name" 


"lim 


<?lim 


''lim 


'lira 


Z[im 


Jlim 


Ks.lira 


3.6yumi,„ 


4.5^m,„^ 


5.8^mii„, 


8.0^m,i„, 




[magAB] 


[magAB] 


[magAB] 


[magAB] 


[magAB] 


[magAB] 


[magAB] 


[magAB] 


[magAB] 


[magAB] 


[magAB] 


SpARCS-0034 


.h 


25.3" 


24.4" 


24.3" 


23.9'-- 


22.5= 


22.2= 


21.4s 


21.2s 


19.7s 


19.68 


SpARCS-0036 


_ h 


25.1'' 


24.4" 


23.7" 


23.5" 


22.7= 


21.5= 


21.2s 


21.18 


19.98 


19.68 


SpARCS-0035 


. h 


25.4'' 


24.9" 


24.3" 


23.6'-- 


24.1' 


23.4' 


22.8S 


22.38 


20.88 


20.48 


SpARCS-0215 


24.8" 


25.1'' 


24.7" 


24.4" 


23.7" 


22.8= 


22.0= 


21.3s 


21.18 


19.58 


19.48 


SpARCS-1047 


25.5" 


25.7" 


25.0" 


24.7" 


23.8" 


23.1'' 


22.3" 


21.6S 


21.38 


19.78 


19.78 


SpARCS-1051 


25.6" 


25.8" 


25.2" 


25.0" 


24.0" 


23.2" 


22.4" 


21.7s 


21.38 


19.88 


19.78 


SpARCS-1613 


25.5" 


26.0" 


25.4" 


24.7" 


24.0" 


23.1" 


22.7" 


22.7s 


22.68 


21.28 


20.98 


SpARCS-1616 


25.1" 


25.7" 


25.1" 


24.8" 


23.5" 


23.3" 


22.7" 


22.6S 


22.48 


21.28 


20.98 


SpARCS-1634 


25.6" 


26.1" 


25.6" 


25.1" 


24.4" 


23.7" 


23.1" 


23.2S 


23.28 


21.68 


21.38 


SpARCS-1638 


25.4" 


25.9" 


25.4" 


25.1" 


24.2" 


23.4" 


22.8" 


23.0S 


23.18 


21.38 


21.38 



MegaCam, Canada-France-Hawaii Telescope (CFHT) 

IMACS, Magellan Telescope 

MOSAIC-II, Blanco Telescope, Cerro Tololo Inter-American Observatory (CTIO) 

WIRCam, Canada-France-Hawaii Telescope (CFHT) 

ISPI, Blanco Telescope, Cerro Tololo Inter-American Observatory (CTIO) 

HAWK-I, Very Large Telescope (VLT) UT4 

IRAC, Spitzer Space Telescope 

Note that for three of the clusters in the Southern hemisphere we do not have li-band data. 



on the SExtractor parameter FLUX_RADIUS measured in the 
Ks-band image, such that q = 0.85FLUX_RADIUS. The factor 
of 0.85 is chosen to optimise the S/N of a source with a circu- 
lar Gaussian PSF-profile. Further we make sure that q is chosen 
such that q > gpsF in all filters. 

This method is applied to measure fluxes in the u - K^- 
bands, but since the IRAC data suffer from a much larger PSF 
we work in a two stage process to incorporate the IRAC fluxes in 
a way that reduces the problems from confusion. We construct 
a 2-stage multi-colour catalogue where we multiply the IRAC 
flux measured in the bigger aperture with the fractional differ- 
ence of the Kj-band flux measured in the small and bigger aper- 
ture. This way we effectively correct the IRAC flux for blend- 
ing with nearby objects by assuming these neighbours have the 
same (Kj-IRAC) colour as the source. For contaminating galax- 
ies this is often the case. To verify that any residual blending in 
the IRAC bands does not affect our results, we repeated the anal- 
ysis while excluding the IRAC data in all SED fits. We find no 
bias in the stellar mass estimates, and even for the lowest masses 
(]VI* = 10'°]VIo), 68% of the estimated stellar masses differ by 
less than 0.05 dex from our fiducial analysis. 

We calibrate the photometric zeropoints on a cata logue basis 
by ma king use of the universalit y of the stellar locus (iHigh et al.l 
l2009l) . We use stellar data from ICovev et al.l (l2007l) . containing 
600,000 point sources selected from the SDSS and 2MASS sur- 
veys. By applying linear colour terms we compare these colours 
to stars measured with the filter sets in the telescope we used. 
Note that these data are especially favourable to calibrate the 
zeropoints using the stellar locus since the amount of galactic 
dust is very low in these fields. We adapt the zeropoints of the 
ugrizSK^ filters to bring the colours of stars in our data in line 
with the reference catalogue. Corrections are typically on the 
order of 0.05 magnitudes. To account for uncertainties in the ab- 
solute zeropoint of IRAC, we included a 10% systematic error 
to the IRAC fluxes. 

After gaussianisation, the background noise in the images 
is correlated between pixels. Therefore we estimate the errors 



on the flux measurements in the stacks of each filter by mea- 
suring the fluctuations in the flux values measured from aper- 
tures that are randomly placed on the images. We take account 
of the non-uniform exposure time over the image stacks. Table 
lA.ll shows an overview of the median 5-cr flux measurements for 
point sources in each filter and each cluster 



Appendix B: Field SIVIF measurements from 
GCLASS 

Thanks to the relatively wide areas that were observed to obtain 
the GCLASS multi-colour catalogues (15' x 15' centred on the 
clusters in the Northern sky, and 10' x 10' for the clusters in 
the Southern sky), these data can also be used to study galaxies 
outside the clusters and hence to measure the SMF of the general 
field. In this appendix we measure the field SMF from GCLASS 
in the redshift range 0.85 < z < 1.20 and compare this to the 
field SMF measured from Ultra VISTA. 

Since the UltraVISTA sample is based on a relatively deep 
(compared to GCLASS) 30-band photometric catalogue, it is 
complete in the mass range (M* > 10'" M©) at this redshift 
range. A comparison between GCLASS and UltraVISTA may 
reveal possible systematic differences in the stellar mass cata- 
logues, and any residual incompleteness in GCLASS. 

To minimise the contamination by cluster galaxies in the 
sample, we use a conservative selection of field galaxies in 
GCLASS. A galaxy is considered as part of the field when it 
is separated from the cluster centre by more than the angular 
distance that corresponds to 1 .5 Mpc at the redshift of the clus- 
ter. Furthermore we require a field galaxy to have a photomet- 
ric redshift |zphot - Zciustei-I > 0.05. After taking account of the 
areas masked by bright stars, this results in a total probed vol- 
ume of the field that is ~ 6 times smaller in GCLASS com- 
pared to UltraVISTA. Since the 10 GCLASS pointings have dif- 
ferent depths, we have to take account of the estimated mass- 
completeness of the detection bands. This is measured similarly 
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Fig.B.l. The UltraVISTA (magenta) versus GCLASS (black) field measurements. The left panel shows the total galaxy population 
in both fields. The middle panel shows the SMF for the subset of star-forming galaxies. The right panel shows the subset of 
quiescent galaxies. Error bars show the 68% confidence regions for Poisson error bars. The grey curves show the 10 contributions 
to the field SMF around the GCLASS clusters, which differ because of cosmic variance due to the small volumes probed in these 
individual fields. Also the fields contribute only down to a particular mass respecting the varying depths of the GCLASS fields. In 
the bottom panels the fractional differences between the two field measurements are shown, given by '^'^^^faavisTA'^^^ ' together 
with the estimated errors. 



as Sect. 13.21 but using a redshift limit of 1.20 in each field (in- 
stead of the individual cluster redshifts). This way we correct 
for Malmquist bias in a similar way as in the l/Vmax weighting 
method. 

Fig. lB.ll shows a comparison of the field SMF measured in 
the GCLASS (black) and UltraVISTA (magenta) surveys. The 
curves are normalised with respect to the total volume subtended 
by these surveys. The grey curves show the contributions to the 
field SMF of the 10 individual GCLASS fields. These contri- 
butions differ between the pointings because their depths are 
different, and also the area surrounding the cluster that is part 
of the field differs. The differences in the grey curves are fur- 
ther caused by cosmic variance. The field in the SpARCS-1047 
image for example is significantly overdense in the redshift bin 
0.85 < z < 1.20. Note however that, when these 10 indepen- 
dent sight-lines of GCLASS are com bined, the uncertainty d ue 
to cosmic variance is greatly reduced (ISomerville et al.ll2004h . 

There is generally a good agreement between the field SMF 
measurements from GCLASS and UltraVISTA, especially at 
the high-mass end. This indicates that there are no substantial 
systematic differences between the two catalogues this study is 
based on. At the low-mass end of the SMF there are some sys- 
tematic differences in both the star-forming and quiescent pop- 
ulation, increasing to several ~ 10% in the lowest mass bins. In 
Sect. 14. II we explained that we corrected the GCLASS cluster 
SMF data by these completeness correction factors. That way 
we can not only compare the cluster and field qualitatively, but 
have a more realistic view on the absolute Schechter parame- 
ters. Note that this additional completeness correction does not 



change any of the qualitative statements in this paper, nor affects 
the conclusions of this paper in any way. 
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